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We investigate possible sky survey modes with the Middle Sized Telescopes (MST, aimed at covering the energy range from 
~100 GeV to 10 TeV) subsystem of the Cherenkov Telescope Array (CTA). We use the standard CTA tools, CORSIKA and 
sim_telarray, to simulate the development of gamma-ray showers, proton background and the telescope response. We perform 
simulations for the H.E.S.S.-site in Namibia, which is one of the candidate sites for the CTA experiment. We study two previously 
considered modes, parallel and divergent, and we propose a new, convergent mode with telescopes tilted toward the array center. 
For each mode we provide performance parameters crucial for choosing the most efficient survey strategy. For the non-parallel 
modes we study the dependence on the telescope offset angle. We show that use of both the divergent and convergent modes results 
in potential advantages in comparison with use of the parallel mode. The fastest source detection can be achieved in the divergent 
mode with larger offset angles (~ 6° from the Field of View centre for the outermost telescopes), for which the time needed to 
perform a scan at a given sensitivity level is shorter by a factor of ~23 than for the parallel mode. We note, however, the direction 
and energy reconstruction accuracy for the divergent mode is even by a factor of ~ 2 worse than for other modes. Furthermore, we 
find that at high energies and for observation directions close to the center of the array field of view, the best performance parameters 
are achieved with the convergent mode, which favors this mode for deep observations of sources with hard energy spectra. 
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1. Introduction 

Imaging Air Cherenkov Telescopes (lACT) detect gamma 
rays using the Cherenkov images of their electromagnetic 
showers developing in the atmosphere. The lACT technique 
has rapidly progressed over the last 20 years (see, e.g., a re¬ 
view i nj ijl) and, with the current generation of lACT instru¬ 
ments 01, it is now the most accurate and sensitive detec¬ 
tion technique in the very high energy gamma-ray astronomy. 
The Cherenkov Telescope Array (CTA), the next generation of 
lACT detectors currently in final stages of design, is expected to 
improve the sensitivity of present observatories by an order of 
magnitude, covering the energy range from a few tens of GeV 
to hundreds of TeV 10 . The experiment will consist of two ar¬ 
rays, one in the northern and one in the southern hemisphere, 
each including a large number of telescopes (~ 30 for north¬ 
ern and ~ 100 southern hemisphere), with a large field of view 
(FOV) of 5°-10° (depending on the telescope type). This will 
allow for several schemes of observation; 

• deep observations - all telescopes pointed onto one object 
(intensive data taking); 
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• normal observations and monitoring - a few telescopes 
oriented towards each of several potentially interesting 
sources; 

• sky scans - all telescopes scan a large area of sky in the 
long-term observations to detect new or transient sources. 

In this work we thoroughly investigate the last one, i.e. sky 
scans. 

On the experimental ground, scans were first performed by 
the HEGRA system of lACTs in search for a TeV gamma-ray 
signal from one quarter of the Galactic plane 0 Ql . At the 
same time, the Whipple telescope collected data from nightly 
calibration scans which covered a sky region 12.5°-wide in dec¬ 
lination 0. Neither of these observations detected TeV emis¬ 
sion. Nevertheless upper limits were derived for a large number 
of sources and the sky-scan techniques and data analysis meth¬ 
ods were developed. Later, during the 1400 h-long survey of 
the Galactic plane conducted by the H.E.S.S. telescope system 
0|T^[Tlt], over a dozen new sources were detected |fl2[]. 

For CTA, an improved Galactic plane survey should be a 
major objective and it will also be capable of performing an 
all-sky survey in unprecedentedly short time at high sensitiv¬ 
ity; the scientific rationale and feasibility of both survey types 
are thoroughly discussed in Q. As also discussed in jT^, 
such surveys can be performed in various modes of observa¬ 
tion, in particular, large number of high-performance lACTs 
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allows for using non-parallel modes with an enlarged FOV. The 
proper adaptation of such a mode for a specific telescope array 
can be a non-trivial task. The optimization of the pointing strat¬ 
egy, taking into account numerous characteristics of an array, 
e.g. distance between telescopes, FOV, energy threshold etc, 
can signihcantly reduce the observation time needed to achieve 
a given sensitivity. 

In this work we consider the array of middle sized telescopes 
(MST) working in various, parallel and non-parallel, modes. 
By performing high-statistics Monte Carlo (MC) simulations 
of the sky-survey observations, we derive for each mode the 
basic performance parameters at both trigger and analysis lev¬ 
els, which then allow us to compare efficiencies of the modes. 
Our study is a part of an intensive work within the CTA Monte 
Carlo Work Package aimed at optimizing the CTA observation 
scheme. Whereas we consider in detail different modes with the 
MST array, independent investigations are currently performed 
for the divergent mode of Large Sized Telescopes (LST) sub¬ 
array and the full CTA array working in divergent modes. 


2. Sky survey modes 


Fig. [T] illustrates possible modes for a large telescope array 
used for sky surveys. The parallel and divergent configura¬ 
tions were considered before in ifisll : below we introduce also 
a novel, convergent mode (note the difference between our ter¬ 
minology and that of OOll . were the parallel mode is referred to 
as convergent). 

The performance of a telescope system operating in the sky 
survey mode depends on the FOV of the system and the time of 
observation needed to achieve a given signihcance level, i.e. its 
sensitivity. 

In the simplest approach, sky surveys may be performed with 
telescopes pointed parallely into the same direction of the sky 
(Fig. m, however, in such a case the FOV of the telescope 
system is highly limited by the FOVs of individual telescopes. 
The FOV of a telescope array can be significantly enlarged by 
slightly deviating the pointing direction of each telescope. In 
the divergent mode, telescopes are inclined into the outward di¬ 
rection, see Fig. [Ib, by an angle increasing with the telescope 
distance from the array center. As explained below, a perfor¬ 
mance improvement for such a conhguration can be expected 
primarily at high energies of primary photons. 

For the divergent configuration, images of gamma rays im¬ 
pinging close the array centre are shifted toward the camera 
edge, which leads to a leakag4i] or complete loss of an event. 
While the larger loss of events is mostly pronounced for the 
lower-energy gamma rays, the leakage effect concerns mainly 
events with higher energies. As a result even if an event is reg¬ 
istered it is poorly reconstructed. On the other hand, orienta¬ 
tion of telescopes in the divergent mode is suitable for efficient 
detection of events with large impact parameter and/or arriv¬ 
ing from directions further from the FOV center (in both cases 
mainly with high energies). 


a) 



b) 



c) 



Figure 1: Three modes of configuration of the telescope system used in the 
sky-suiwey scans: a) normal (parallel) mode; b) divergent mode; c) convergent 
mode. 


Qualitatively, one can expect that those negative effects can 
be reduced for the opposite orientation, i.e. with outer tele¬ 
scopes inclined toward the array center, see Fig. A quan¬ 
titative comparison of the performance of the three modes and 
a related issue, i.e. an optimal value of the offset angle (giv¬ 
ing the amount of the difference of the pointing directions, as 
dehned below), appears crucial for planning the most efficient 
survey strategy. 

3. MC simulations 

For all three modes, we simulate the response of the tele¬ 
scope array to the Extensive Air Showers (EAS) induced by 
gamma rays and proton background. To simulate the develop¬ 
ment of EAS we use CORSIKA 6.99 code lITillTsIl . used as 
a standard in CTA. We simulated 2.1 x 10’ gamma rays and 
3.8 X 10* proton event^ - both with energies between 30 GeV 
and 10 TeV generated from differential spectra with the spectral 
index F = -2.0. However, in our analysis, we use event weights 
corresponding to spectra with F = -2.57 for gamma rays and 
F = -2.73 for protons. Gamma rays are simulated from a point¬ 
like test source with the direction defined by the Zenith angle 
Za=20° and the Azimuth Az=180° measured with respect to 


* The effect of cutting off an image at the camera edge. 


^including the number of re-used showers 
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Figure 2: Geometric layout of telescopes used in our simulations. The circle 
shows the simulated area defined by the maximum impact pai'ameter for gamma 
rays. 


the magnetic North. The background proton showers are simu¬ 
lated isotropically with directions within a 10° half-angle cone 
(larger than the FOV of all considered modes) centered on the 
direction of the gamma-ray source. We set the maximum im¬ 
pact parameter for gamma rays to 1000 m and for protons to 
1500 m. The detector array is assumed to be located at the 
Namibian (H.E.S.S.) site at the altitude of 1800 m a.s.l. 

The response of the telescope array is simulated with the 
CTA sim_telarray code ifla fl^ . We use the MST sub¬ 
array of the CTA array E from the so-called production-1-, the 
subarray includes 23 telescopes with positions shown in Eig.|3 
The direction of the central telescope no. 5 is always approxi¬ 
mately in the center of the FOV of the array (a slight displace¬ 
ment may occur due to the presence of telescopes no. 12 and 15, 
which break the symmetry); then, this direction is used to define 
various configurations and we refer to it as the axis of the tele¬ 
scope system. The non-parallel modes are defined by altering 
the pointing-directions of all remaining telescopes with respect 
to this axis (outwards for divergent mode, inwards for conver¬ 
gent mode) by an offset angle, a, givenbj0a = k-t/ter^tei/2rniax, 
where k is the scaling parameter (=1 or 2, see below), t/tei = 8° 
is the FOV diameter of a single telescope, rtei is the distance 
of the telescope from the array center and r^ax = 480 m is 
the largest distance (i.e. the one of telescopes no. 20-23). For 
both non-parallel modes we consider two scales of the angular 
spread between telescopes, with k - I and k - 2, referred to as 
single and double scaled modes, respectively. Then, in total we 
consider five configurations: the normal mode (N), the single- 
scaled divergent mode (D), the double-scaled divergent mode 
(2D), the single-scaled convergent mode (C) and the double- 
scaled convergent mode (2C). For example, the telescopes no. 6 
(with rtei = 120 m) and no. 9 (with rtei = 170 m) are offset by 


^except for the most distant telescopes 20-23, for which we use a = 3k ■ 

‘^tel ■ r tel/8rmax- 



Figure 3: Telescope pointing directions shown in the Zenith (radial) and Az¬ 
imuth (transversal) polar coordinates for mode D (a), C (b), 2D (c) and 2C (d). 
For each mode, telescopes no. 20 to 23 are indicated with different markers 
for better illustration of the difference between the divergent and convergent 
modes. 


Sim telarray input 

Input value 

Telescope type 

MST (Type 2) 

Dish diameter 

d = 12 m 

Focal length/diameter 

II 

Camera Field of view 

FoV = 8° 

Pixel size 

0.18° 

Photomultipliers 

bi-alkali 

quantum efficiency 

QEpeak = 25.7% 

Telescope trigger 

Min. 4 pe in each of 

threshold level 

3 neighboring pixels 

Min. trigger multiplicity 

2 telescopes 


Table 1: Parameters assumed in sim.telarray for our simulations. 


a - 1° and 1.4°, respectively, for k - I and by 2° and 2.8° for 
k — 2. The orientation of telescopes in non-parallel modes is 
shown in polar coordinates in Fig. [3] 

We emphasize that the term ’’convergent” used throughout 
this work refers to pointing scheme where the individual tele¬ 
scope axes cross in an arbitrary point of a shower core. Note 
that in the C mode the telescope axes converge approximately 
at a distance of about 7 km, close to the shower maximum for 
TeV showers while in the 2C mode the point of convergence is 
well below the maximum of all showers. 

The basic technical parameters of the telescopes used in our 
simulations follow the standard settings of the production-1 (the 
hrst CTA MC mass production; see Fig. 18 and Chapt. 8 of ll^ 
or Chapt. 6 of 11^ 1 and are given in Table [1] We note that 
the currently developed MC mass production prod-2 introduces 
several technical modifications (e.g. longer read-out intervals 
and trace integration following the time gradients) which may 
improve the efficiency of detection of events with impact pa- 
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Figure 4: The number of registered gamma-ray events as a function of the impact parameter measured from the array center to the shower core. The top panels are 
for E < 70 GeV (the array threshold) and the bottom panels are for E > 1 TeV. For both energy ranges, 0ofF = 0° (on-axis), 4.9° and 6.3° from left to right. In 
each panel the modes are represented by: N - solid lines with filled circles, D - dotted lines with filled triangles, 2D - dotted lines with empty triangles, C - dashed 
lines with filled squares and 2C - dashed lines with empty squai'es; the same styles of markers and lines are used for the modes throughout all figures in this paper. 
In all cases the same number of simulated events were used, so the diflFerences in the number of entries (y-axis) between the modes in each panel correspond to the 
expected difference of trigger rates. 


rameter above ~ 400 m. 

Obviously, sky scanning surveys are aimed at the detection 
of unexpected and non-resolved sources, which may be located 
anywhere in the scan area. Then, we take into account source 
locations in various regions of the FOV of the array. To this 
end we simulate detection of source with different angular dis¬ 
tances, (pos - f(Za, Az), from the array axis by moving the FOV 
of an array with respect to the simulated source direction. As¬ 
suming the rotational symmetry of an array FOV in (Za, Az) 
plane, we set Za to zenith angle of the array axis direction and 
we change only the Az component. Specifically, we consider 
test sources with different 0off. with 0oft between 0.0° and 7.0° 
and we use the step of A^off = 0.7°. The case of the source with 
(poS - 0.0° is referred to as the on-axis source and all cases with 
(pos + 0.0° as off-axis. 

Our analysis is based on Hillas image parameters l lSll ob¬ 
tained from Sim_telarray simulations using the read_cta 
code, being a part of the CTA simulation package. Images are 
cleaned using the tail-cuts of 5.5/11 photo-electrons (pe). 

The shower direction and energy reconstruction used in this 
work follows the CTA standard baseline methods as described 
in Sect. 3 of iflTll . The shower direction is determined with clas¬ 
sical stereo reconstruction which uses Hillas image parameters 


obtained from individual telescopes and then transformed into 
the common reference frame (CRF) of coordinates. In our ap¬ 
proach CRF is defined by the plane perpendicular to the array 
axis. The use of CRF is a general method yielding proper re¬ 
construction for both the parallel and the non-parallel pointing 
schemes. The shower parameters, which result from a stereo re¬ 
construction carried out in CRF, then are used in the energy re¬ 
construction. The energy reconstruction, in turn, uses the corre¬ 
lation matrices between energy, image SIZE and reconstructed 
impact parameter. Parameters of correlation matrices are stored 
in the so-called lookup tables (LUTs). In this work, to avoid the 
possible biasing of the results, the energy reconstruction is done 
in two steps using two separate event samples. First, with the 
gamma-ray training sample, we create LUTs. Then, the LUTs 
are used to determine the energy of both gamma rays and pro¬ 
tons in the independent test samples. As the stereo parameters 
stored in LUTs depend on the observing mode, we use LUTs 
specific for each studied mode. 


4. Results 

We compare the basic parameters used to describe the per¬ 
formance of ground-based gamma ray detectors, cf. Hi], sep- 
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arately at the trigger and analysis levels, i.e. before and af¬ 
ter gamma/hadron separation cuts. The analysis level param¬ 
eters allow for a straightforward comparison of the overall per¬ 
formance of the modes, however, they depend on the chosen 
gamma/hadron separation method and then they can be im¬ 
proved if optimized procedures are developed. In turn, the 
trigger level parameters provide information independent of the 
chosen analysis method and then they are directly related to the 
detector physics. 

4.1. Trigger level parameters 

4.1.1. Impact parameter distributions 

We first illustrate basic effects underlying the difference of 
performance discussed in following sections. The optimization 
of the performance should maximize the number of gamma-ray 
events which (1) are detected, and then (2) have properly recon¬ 
structed arrival direction and energy. While (1), i.e. efficient de¬ 
tection, is mainly determined by the single telescope response, 
(2) improves with a larger number of triggered telescopes (re¬ 
gardless of the analysis method). In Figs.|4]and|5]we show the 
relevant distributions of the number of registered events, 
and the number of triggered telescopes, Nt^u as functions of the 
impact parameter, IP, measured from the array center to the 
shower core. 

Low energy events with energies around the energy threshold 
of the MST array, E ~ 70 GeV, are typically registered by 
a single telescope at distances not exceeding the radius of the 
Cherenkov light pool on the ground, ~ 120 m, which then for a 
whole telescope array corresponds to the IP ~ 200 m. As seen 
in Figs.lUa-c), for a large cpoff > 5° (i.e. larger than the radius of 
FOV of a single MST) the use of non-parallel modes has a huge 
effect in the low energy range; all non-parallel modes give n^ 
by 2 orders of magnitude larger than the N mode. Furthermore, 
in this regime (i.e. with low E and large (poff), the divergent 
modes are notably better than the convergent modes especially 
at larger IP. For sources located on-axis, ntr(IP) is almost the 
same for the N, D and C modes, whereas for the 2D and 2C 
modes «tr is much smaller, which shows that in this regime (low 
E and small 0off) the leakage effect is important only for larger 
telescope offset angles (a). 

High energy events, in the TeV range, produce a large density 
of Cherenkov photons, which is sufficient to trigger a telescope 
at distances significantly exceeding the light-pool radius, hence 
on average IP > rmax- At high energies, both divergent modes 
give much larger n^ than the C, 2C or N mode regardless of <poff, 
except for small 0off for which N and D give similar (largest) 
values of n^, see Figs|3d-f). 

As also seen in Fig.|4] at all energies the single scaled modes 
give much larger Wtr than the corresponding double scaled 
modes up to (poS ~ 5°. However, at large (pos > 5° the dou¬ 
ble modes are more efficient and especially at high energies the 
2D mode gives a much better detection efficiency than other 
modes. 

The number of triggered telescopes is much smaller for low 
energy events, on average Atei ~ 2-3 for all observation modes 
and any location of the source, than for high energy events. 






Figure 5: The average number of telescopes triggered by gamma-ray events as 
a function of the impact parameter (measured from the array center) for the five 
modes, represented by the same line and marker styles as in Fig.|4] Top panels: 
E <10 GeV, bottom panels: E >\ TeV, left: ^off = 0°, right: (^off = 4.9°. 


However, for the latter it strongly depends on both the source 
location and the mode. Namely, for an on-axis source (Fig.|5};) 
and for small IP (< 300 m), Atei is much smaller for modes D 
and 2D (for the latter by over a factor of 2) than in other modes, 
which is an obvious consequence of the loss of the event in the 
outer, most diverged telescopes. A more complex effect oc¬ 
curs for off-axis sources, for which in some configurations the 
displacement of the image from the camera center due to the 
increasing IP balances the (opposite) displacement due to the 
increasing 0off in some ranges of IP, however, it results in a sig¬ 
nificant increase of Atei only in mode C (see Fig. Si). 

Briefly summarizing the above, we can expect the largest 
trigger rates in mode D (for smaller (pos) or 2D (for larger (pos) 
and the most precise reconstruction in mode C. The overall per¬ 
formance results from the combination of these two features 
and then the above basic properties do not allow to assess which 
mode can be most suitable for a specific type of observation. To 
this end a extensive comparison of the standard lACT perfor¬ 
mance parameters is needed with use of a possibly most effec¬ 
tive analysis methods. 

4.1.2. Collection areas and trigger rates 

The collection area of an lACT array, A{E), can be ob¬ 
tained from MC simulations from the number of triggered and 
simulated events, n^ and Wsim, and the simulated area, S, as 
A{E) - [ntr(F)/«sim(F)] X S. Fig. |6] presents the E- and (^off- 
dependence of A. Essentially it reflects the dependence of n^ 
on E and (pos discussed above. In particular, 

(i) the largest values of A correspond to (pos = 0°; then, A 
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Figure 6: Trigger-level collection area as a function of E (top panels) and 
(bottom panels) for the five modes, represented by the same line and marker 
styles as in Fig.|4] In the top panels the black lines and markers are for <j)„g = 0° 
and the grey ones for 0off = 4.9°. Panel (c) is for E = 10 GeV and panel (d) for 
E = 5 TeV. 


decreases with increasing (pos, but the shape of Aipos) signifi¬ 
cantly differs between the modes. The most rapid decrease oc¬ 
curs in mode N, the flattest A(0off) correspond to double scaled 
modes, for which A is larger than for single scaled modes at 
> (4 - 5)°. 

(ii) for small 0off and low E (~ 70 GeV), A in double scaled 
modes is roughly half that of single-scaled modes. 

(iii) for small (pos and high E 5 TeV), A is approx, by a factor 
of 2 larger in divergent modes than in convergent modes. 

The collection area provides information independent of the 
energy spectrum of a source. To illustrate efficiency of the 
modes in a realistic data-taking, we use as a reference the spec¬ 
trum measured from the Crab nebula around 1 TeV, Fy{E) - 
2.79 X 10-"(£/lTeV)-2-5^ cm-^s-'TeV-* (see Chapt. 8 in 11), 
and we calculate the total trigger rates, Ry, by integrating over 
E the differential trigger rates, Ry(E) - Fy(E) x A{E). The re¬ 
sulting Ryiipos) is shown in Fig.|7] The major contribution to Ry 
comes from low E events, then the dependence on cpoff is close 
to that of A((poff) at low E (shown in Fig.|6};). 

We also find the total trigger rate for protons, Rp, by 
integrating over E the differential trigger rate, Rp{E) = 
f dQ(a))Fp(E) X Ap(E,a>), where Fp(E) - 8.72 x 

10“®(£'/lTeV)“^'^^^ cm“^s“'sr''TeV“' (i.e., at 1 TeV compara¬ 
ble with proton flux obtained with BESS-TeV experiment 1201] ), 
Ap is the collection area for protons and Q is the simulated 
proton cone angle related with the half opening angle, w, by 
Q(w) = 27r(l - cos w). 

Comparing the values of Rp given in Tablel^with Ry in Fig.|7] 



Figure 7: The total trigger rates for gamma rays, Ry, for a Crab-like source (see 
text) as a function of (pog for the five modes, represented by the same line and 
marker styles as in Fig.|4] 


Total background proton rates [protons/s] 

N 

D 

C 

2D 

2C 

7888 

7446 

8224 

6633 

8032 


Table 2: Total background proton rates, Rp, for the studied modes. 

we see that for an on-axis source Rp are larger by three orders of 
magnitude even for the most efficient (in gamma events) modes 
N, D and C, and the ratio of Ry to Rp decreases with increasing 
(poS- Notably, Rp in the convergent modes is by ~ 15 % larger 
than in the divergent mode^. 

We emphasize that in our background simulations we do 
not take into account other components, e.g. helium, electrons 
etc., and thus some performance parameters, e.g. absolute val¬ 
ues of sensitivities, could differ from more realistic estimations 
for CTA. On the other hand, it should not affect significantly 
the relative comparison of different modes. We briefly discuss 
expected changes in sensitivities due to contribution of other 
background components at the end of Sect.|5] 

4.1.3. Acceptance 

As noted above, the A(0off) dependence is different in the 
considered modes, in particular, it is much flatter in modes 
2C and 2D, whereas in modes N, D and C it is strongly en¬ 
hanced at (poff <5°. To compare the detection efficiency in 
the total observed area, we define the acceptance, Acc(£’) = 
f A(E, (poff)dQ, where dQ = 2n sin (p^s x d^off (assuming the 
variation of A{E, (pos) across the FOV depends exclusively on 
(pos)- As df2 increases with the offset angle (pos, the relatively 
stronger contribution to Acc{E) comes from the regions with 
larger (pog. In other words Acc(£’) can be used as a measure 
of the performance of the array for an isotropic population of 
sources within the FOV 

As we see in Fig. [8j the difference between Acc in different 
modes is larger at high energies, where also larger Acc corre¬ 
spond to the divergent modes and the largest (for mode 2D) and 


''note that apart from a potential advantage in a better significance, arrays 
triggering less background events are preferred for easier acquisition of data 
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Figure 8: Acceptance at the trigger level as a function of energy for the five 
modes, represented by the same line and marker styles as in Fig.|5] 


smallest (for mode C) acceptances differ by over a factor of 2. 
At low energies the differences are smaller, but also here the 
largest Acc corresponds to modes D and 2D and the smallest to 
mode 2C. 

4.1.4. Reconstruction 

In this section we compare the accuracy of the direction and 
energy reconstruction in different modes. The geometric re¬ 
construction method used for lACT relies on the Hillas image 
parameters, describing geometrical properties of an image in 
a telescope camera plane. The arrival direction of a primary 
particle is estimated from the intersection of the major axes of 
Hillas ellipses of images from each pair of triggered telescopes 
(averaged over all pairs if more than two telescopes triggered). 
The quality of the direction reconstruction is given by the an¬ 
gular resolution, defined by the containment radius, ry, of the 
area in the camera plane enclosing a fixed fraction f of the re¬ 
constructed directions (for details see Sect. 4.5 in 12111 '). 

The angular resolution for / = 68% is shown in Fig. |9l The 
difference between the modes results mostly from two effects. 
First, the quality of the direction reconstruction improves with 
the increasing number of telescopes triggered by an event, as a 
larger number of major axis crossing points tends to reduce the 
reconstruction inaccuracies. Secondly, the angular resolution is 
spoiled by the leakage effect (the truncated images typically do 
not contribute to improvement of reconstruction). We note that 
both effects add up and yield a much better angular resolution 
in the convergent than in the divergent modes; in the former Atei 
is typically much larger (see Fig.|5]) and in the latter the leakage 
effect is typically stronger, especially at high energies. An ad¬ 
ditional, but weaker effect involves different shapes of images, 
with more elongated images allowing for a more precise deter¬ 
mination of the image major axis and, then, a better estimation 
of the particle arrival direction. This effect results in a slightly 
better angular resolution in mode C than in D at small E and 
large (^off, for which mode C has a larger average IP (implying 
more elongated images), see Fig.|5] 

In standard analysis of parallel mode configuration the angu¬ 
lar resolution typically improves with increasing E, for which 
Hillas ellipses are better constrained. Such an improvement in- 




Figure 9: Angular resolution before background suppression cuts for the 68 % 
containment radius as a function of E (a) and 0off (b), for the five modes rep¬ 
resented by the same line and marker styles as in Fig.|4| In panel (a) the black 
lines and markers are for = 0° and the grey ones for 0off = 4.9°. In 
panel (b) the black lines and mai'kers are for £ = 3 TeV and the grey ones for 
£ = 70 GeV. The lack of several points in panel (b) is due to an insufficient 
statistics to compute the angular resolution. 


deed occurs in mode N and C for an on-axis source (Fig. |9^). 
However, in other modes the angular resolution worsens with 
increasing E, above ~ 500 GeV, and we note that this is mostly 
due to the leakage effect, by which modes D, 2D and 2C are 
strongly affected at TeV energies. As a result, large differences 
occur in this large energy range, in particular at £ > 3 TeV 
the best resolution of ^ 0.07° in mode C is by over a factor 
of 3 better than the resolutions in modes 2D and 2C and twice 
better than the resolution in mode D. Notably, the resolution in 
modes N, D and C significantly worsens with the increase of 
0off. whereas for 2D and 2C it is less sensitive to 0off and at 
‘/’off ^5° the double scaled modes have a better resolution than 
single scaled ones. At low energies, ~ 100 GeV, all modes have 
a poor resolution, > 0.2°, and it is less dependent on both (pos 
and the mode. 

The estimation of the primary particle energy relies on its 
correlation with the SIZE of an image (the total number of 
photo-electrons in image) and the impact parameter. SIZE is 
obtained by integrating the light collected in each pixel belong¬ 
ing to an image and the impact parameter is estimated with use 
of DIST Hillas parameter (distance between the reconstructed 
source position and the center of gravity of an image), to which 
it is proportional. The energy reconstruction quality of lACT 
telescope systems is conventionally described by two quanti¬ 
ties, the energy resolution and the energy bias. The first one 
is defined by the value of the standard deviation of a Gaussian 
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Figure 10: Top panels show the trigger-level energy resolution as a function of 
E (a) and (b) and bottom panel (c) shows the trigger-level energy bias as 
a function of for the five modes represented by the same line and marker 
styles as in Fig.|4] In panel (a) the black lines and markers are for = 0° and 
the grey ones for = 4.9°. In panels (b) and (c) the black lines and markers 
are for E = 3 TeV and the grey ones for £ = 150 GeV. 


fit of the distribution of (E^ec - E)IE, where E^^c and E are the 
reconstructed and true (i.e. simulated) energy, respectively. 

The bias is the mean of the fitted Gaussian. The differences in 
energy biases - the most apparent at energies close to the array 
threshold - reflect the possible changes of the effective energy 
threshold among the studied modes. The value of the effective 
energy threshold, being usually taken as the first point of the 
reconstructed energy spectrum, is relevant to properly identify 
the spectrum of observed source in the possibly largest energy 
range. 

As seen in Fig. (TO] qualitatively, the energy resolution shows 
the same dependence on the mode, E and 0off as the angular 
resolution and this results from the same effects, involving the 
number of telescopes which are triggered and survive the image 
reconstruction (which improves the geometrical reconstruction) 
and the leakage effect which spoils the energy reconstruction 
especially in non-parallel modes (except for mode C) at high 
energies. Apart from this, the standard effects affect the energy 
dependence, namely, the Hillas ellipses are better defined for 
higher E allowing for a more precise estimation of the center 
of gravity of an image, and hence a more accurate reconstruc¬ 
tion of IP, furthermore, higher energy gamma rays give larger 
SIZES, less sensitive to fluctuations. In all modes the best en¬ 
ergy resolution is obtained at several TeV and the overall best 
value of AEjE ~ 0.05 of the mode C is by a factor of ~ 3 better 
than that in mode 2D. 


Also the energy bias is the lowest for mode C, except for 
large 0off (see Fig. [TObi. In all modes, the bias is negative at 
high E and positive at low E. It has the lowest absolute value, < 
0.05, at high E and low 0off, and increases with the decrease of 
E and increase of 0off- The most biased energy reconstruction 
occurs for D mode. At energies of 150 GeV and at offsets 
of > 4° the expected bias for D mode is by 60% larger in 
comparison with bias for both convergent modes. 


4.2. Analysis level parameters 


4.2.1. Gamma/hadron separation and sensitivity 

To separate the gamma-ray signal from the dominating 
hadronic background we apply the standard gamma/hadron 
separation method with shape and direction cuts. Specifically, 
we use cuts on the mean scaled WIDTH (MSW), mean scaled 
LENGTH (MSL) and 6^ (see 122 ] for definitions), which are the 
most effective separators. To improve the background rejection 
we also dynamically optimize the cut on the minimal number of 
telescopes used for reconstruction of an event. We do not apply 
any nominal distance cuts (in e.g. the DIST parameter). 

Cuts are optimized to provide the best sensitivity. To avoid 
biasing the results with optimization on statistical fluctuations, 
we first tune the cuts on separated MC sub-samples (train sam¬ 
ples) of gamma and background events. Using such optimized 
cuts, we calculate sensitivities on independent MC test samples. 
An optimization of cuts with train samples is performed inde¬ 
pendently for each observational mode, each energy bin and for 
each value of (pog. Note that we calculate the point source sensi¬ 
tivity, however, the term sensitivity is used throughout the paper 
for brevity.. 

Sensitivity is the key parameter describing the performance 
of a telescope array. It gives the magnitude of the flux observed 
from a source in a given time at a particular significance level; 
we assume the usual level defined by the gamma-ray signal ex¬ 
ceeding the background estimate by 5 standard deviations (cr) 
and we find sensitivities for an observation time of 50 h. 

In our analysis we use the definition of differential sensitiv¬ 
ity which follows from the Li & Ma significance formula (see 
eq. (17) in commonly applied in lACT experiments, in¬ 
cluding CTA. We also mply the standard CTA conditions for 
a significant detection i.e. gamma-ray signal exceeding 5% 
of the remaining background and at least 10 events in excess, 
Aex- To calculate the differential sensitivity we divide the re¬ 
constructed energy range into five bins per decade and we in¬ 
tegrate gamma and proton differential rates (see Sect. 14.1.11) 
to get ON-axis counts, Aon, and OFF-axis counts, Aqff, for 
each simulated (pog. In our analysis, the off-axis events are 
binned in 11 rings around the array FOV centre. Specifically, 
rings (angular bins) are defined by 11 consecutive ranges of 
offsets: [0.0° - 0.5°], [0.5° - 1.0°], [1.1° - 1.7°], [1.8° - 2.4°], 
[2.5°-3.1°], ... , [6.7° - 7.3°], i.e. off-axis angular bins are 
centered around each (pgg. We checked that the contribution of 
background events with the real offset > 10.0° and the recon¬ 
structed offset < 7.0° is negligible even for modes 2C and 2D. 

In FigHT] we present the sensitivity as a function of i^off in 
several energy bins. Note that for mode N the number of events 
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Figure 11: Differential sensitivity as a function of i^off in four bins of the reconstructed energy bins, given in the panels, for the five modes represented by the same 
line and marker styles as in Fig.|4] Standard CTA conditions for a significant detection are used (5 it significance in 50 h with gamma excess at a level of at least 5% 
of the remaining background). Sensitivity is given in standard HEGRA Crab Units, IC.U. = 2.79 X cm^^s^’TeV^* . 


registered at (ftos ^ 4° is too small to satisfy the detection con¬ 
ditions specified above, so in this mode we formally obtain an 
infinite required flux outside of the nominal FOV of an individ¬ 
ual MST. 

To increase the statistic of proton events, all the events from 
each ring are selected. The total number of events selected in 
a ring is scaled to account for the difference of area between 
the ring and standard circular OFF regions. In each ring, the 
ratio, T], of the ON-source to the OFF-source proton events used 
in background estimation, is set to 1/5. Then, the detection 
conditions are given by A^ex = A^on - > 10 and S 

5%A^bkg, where A^bkg = 

The sensitivity of a telescope system is mainly determined by 
the effectiveness of gamma/hadron separation cuts and the (trig¬ 
ger) efficiency of gamma-events detection; obviously, the dif¬ 
ferential sensitivity depends also on the accuracy of the energy 
reconstruction. We note that the efficiency of separation cuts at 
given E and 0off differs significantly between the modes, with 
differences of the fraction of rejected proton events exceeding 
an order of magnitude. The efficiency of the 9^ cut is directly re¬ 
lated with the angular resolution and the efficiency of the cut on 
the number of triggered telescopes depends on A^ei (see Fig. |5] 
and Fig. [T2| and, then, these separation cuts are particularly 
efficient in both C and 2C modes. Since for the convergent 


modes A^tei rapidly increases with the energy (see Fig. fT^ . the 
cut on the number of triggered telescopes in convergent modes 
becomes particularly efficient above 300 GeV. For those en¬ 
ergies and (/loff > 4°, the sensitivities of convergent modes are 
comparable or even better than those of the corresponding di¬ 
vergent modes, see Figs. fTTT' c-dl. In mode 2C we additionally 
note highly efficient MSW cuts at large (pos. 

The separation cuts in the divergent modes are typically less 
efficient, nevertheless, the sensitivity of mode 2D is much better 
than those of other modes for large (p^s, while mode D has a 
good sensitivity at small (pos, in both cases due to high trigger 
efficiencies of gamma events. 

Summarizing this section, we emphasize two properties cru¬ 
cial for our final conclusions: 

(i) a high sensitivity of mode C for 0off = 0, resulting from 
highly efficient cuts in both the direction and the number of 
telescopes (see Fig. fTZb . 

(ii) a high sensitivity of mode 2D at large resulting from 
the high trigger efficiency. 

As a final remark we emphasize that the sensitivities reported 
in this section serve only the purpose of comparing the differ¬ 
ent pointing configurations. The more advanced gamma/hadron 
separation methods are currently deployed in CTA on a more 
realistic background profiles. The sensitivities provided with 
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Figure 12: The average number of telescopes triggered by gamma-ray events as 
a function of 0off for the five modes, represented by the same line and marker 
styles as in Fig. llll Panels (a)-(d) correspond to energy ranges used in panels 
(aHd) in Fig.[n] 


those methods give a more adequate expectations of the abso¬ 
lute sensitivities of the CTA instrument. 

4.2.2. Collection areas and acceptance 

Rejection of gamma-ray events in gamma/hadron separation 
reduces the collection areas at the analysis level with respect to 
the trigger level by a factor of ~ 2 for ipog = 0° and by a fac¬ 
tor of several for 0off > 4°, see Figs. IhlandflBl This implies 
a relatively weaker contribution from large 0off to Acc at the 
analysis level and, furthermore, the decrease of Acc by a factor 
of several in all modes, see Figs. Island [T4l However, the 2D 
mode has the largest collection areas at large (poff and also the 
largest Acc in the whole energy range, similarly as at the trigger 
level. We also note that the fraction of rejected events is differ¬ 
ent in various modes at low energies, however, at £ > 1 TeV 
it is similar in all modes and hence at high energies the ratio of 
their collection areas or Acc values is roughly the same at the 
analysis and trigger levels. 

4.2.3. Quality of reconstruction 

The background rejection is accompanied with a large rejec¬ 
tion, exceeding a factor of 2 in magnitude, of gamma-ray events 
(c.f. Fig.l^and Fig.fOll. Major part of excluded shower images 
are of low quality, thus their rejection improves the reconstruc¬ 
tion efficiency. Consequently, the reconstruction performance 
of an array at the analysis level is in principle at least as good 
as before introduction of the background suppression cuts. 

In this section we briefly consider changes in reconstruction 
efficiency of studied modes at the analysis level. We do not 



Figure 13: Analysis-level collection area as a function of E (top panel a) and 
cpoS (bottom panels) for the five modes, represented by the same line and marker 
styles as in Fig.|4] In panel (a) the black lines and markers are for (p^ff = 0° and 
the grey ones for (paB = 4.9°. Panel (b) is for £ = 150 GeV and panel (c) for 
E = 5 TeV. 


discuss differences in energy bias as they are comparable with 
those before introduction of gamma/hadron separation cuts (for 
details see Sect. l4.1.4l l. 

By virtue of used cuts the large leakage effect noticed in 
Sect. 14.1.41 which mainly affected divergent arrays, is signih- 
cantly reduced. The resulting improvement in angular resolu¬ 
tion is the most noticeable in the high energy range > 2 TeV 
for 2D mode where for foS = 4.9° it is comparable with N and 
it is even better than 2C (see Fig. ITSk). Interestingly for both 
convergent modes the effect of leakage being the lowest before 
use of analysis cuts, is not further reduced at the analysis level 
- the angular resolutions corresponding to C and 2C do not sig¬ 
nificantly change. 

Similar improvement is noticed in energy resolution of arrays 
set in divergent modes, see Fig.[T5](b and d) (c.f. with Fig.fTOlla 
and b)). For example, at hxed energy of 3 TeV, the energy reso¬ 
lutions for D and 2D modes are improved by almost 25 %. For 
other modes, the energy reconstruction after gamma-selection 
cuts is also more accurate. However the decrease of the corre¬ 
sponding values of energy resolution is less than 10 %. Note 
that the C mode which provides the best energy resolution be¬ 
fore gamma-selection cuts (see Fig. [10] (a and b)) remains the 
best also after use of cuts. 

In this section we present results of sky scans performed with 
the five modes on a hxed sky area of 20° x 12°, see Fig.[T6| For 
modes 2D and 2C we use two times fewer pointings, taking 
into account the large effective FOV and the signihcantly better 
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Figure 14: Acceptance at the analysis level as a function of simulated energy 
for the five modes, represented by the same line and marker styles as in Fig.|4] 




Figure 15: Analysis-level angular and energy resolutions, for the five modes 
represented by the same line and marker styles as in Fig. |4] Panel (a) shows 
the angular resolution for the 68 % containment radius as a function of the 
simulated energy. Bottom panels show the energy resolution as a function of 
the simulated energy (a) and (b). The black lines and markers in panels 
(a,b) are for ipoS = 0° and the grey ones for ipoS = 4.9°. In panel (c) the black 
lines and markers are for E = 3 TeV and the grey ones for £ = 150 GeV. 


sensitivity at large (f>oS in these modes. The total exposure/scan 
time is the same for each mode. To roughly compare the scan 
efficiency of each studied mode, we calculate the mean sensi¬ 
tivity, S mode, achieved during the scan time of 150 h. 

Technically, to evaluate S mode we compute the total gamma- 
ray exposure and the total background rate for each sky coor¬ 
dinate from all the sky pointing positions. Then for each co¬ 
ordinate we calculate the sensitivity (flux fraction) satisfying 
the conditions described in Sec. 14.2.11 S mode is Anally com¬ 
puted by averaging the achieved sensitivities over the whole 


map of coordinates. To avoid possible large contributions from 
regions where the flux fractions are very large (to achieve the 
significance of 5 cr), the averaging is carried out with a use 
of the harmonic mean. The relative statistical uncertainty, 
A5 mode/5 mode, of S mode, estimated with an RMS of the dis¬ 
tribution of significance over a whole sky map, is of order of 
0 . 2 . 

4.2.4. Sky scans 

In all modes, we And the best mean sensitivities in the energy 
bin between 300 and 480 GeV and Fig. [T6]presents the sensi¬ 
tivity map for this bin; the mean values of S at higher and lower 
E are given in Table [3] 


Mean scan sensitivities [% C.U.] 

Energy [GeV] 

N 

D 

C 

2D 

2C 

120-190 

2.5 

2.0 

3.5 

2.0 

3.6 

300-480 

2.1 

1.6 

2.3 

1.4 

1.9 

1900-3000 

4.2 

2.3 

3.0 

2.3 

3.2 


Table 3: Mean sensitivities in sky scans performed with the five modes for the 
pointing layouts presented in Fig. EH 

We see that at £ = (300 - 480) GeV the best scan results are 
obtained with the 2D mode. Its 820 is by 33% better than 
S N, which implies that a scan can be performed with mode 2D 
in a time shorter by a factor (S 2 ^/ 820 )^ ~ 2.3 than with mode 
N. Also for modes D and 2C, 5 d = 1.6% C.U. and S 2 C - 
1.9% C.U. indicate reduced sky-scan times, by a factor of » 
1.7 and ^1.2 with respect to mode N. The worst mean scan 
sensitivity corresponds to mode C. 

Properties similar to the above occur also at low and high E, 
although here 52 d - So- However, we regard the comparison 
of mean 8 at the intermediate E as the most informative. 

We emphasize that the mean scan sensitivities shown above 
resulting from a simple analysis methods and are used only 
for relative comparisons of modes. The more advanced analy¬ 
sis techniques with introduction of a more realistic background 
profiles could lead to a different absolute values of sensitivities. 
On the other hand, the possible contribution of the other back¬ 
ground components is expected to not too significantly change 
the relative differences between mean sky sensitivities obtained 
by the modes. By virtue of a very approximate method of using 
gamma rays with different offsets from the camera center as a 
proxy for diffuse electrons (being one of the most irreducible 
background components) we estimated the possible changes of 
sensitivities due to electrons contribution. The resulting de¬ 
crease of relative differences between the mean sky sensitivi¬ 
ties does not exceed ~ 10% and consequently the final results 
presented in this work, at least qualitatively, hold. 

5. Summary and discussion 

Using the CTA MST subarray, we studied various modes of 
observation with a large array of lACTs. Our results should 
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Figure 16: Sensitivity at energies between 300 and 480 GeV in sky scans performed with the five modes. The sky area of 20° x 12° is covered with the same 
exposure time of 150 h for each mode. We use a simple multiple-row scan strategy to uniformly expose the investigated sky region. Because of the relatively larger 
effective FOV of modes 2D and 2C, the number of pointings for these modes is reduced by a factor of 2. The mean scan sensitivity achieved after the total exposure 
time of 150 h is given for each mode for a robust compaiison of sky-survey efficiencies. 


help in the development of an optimal observing strategy, mak¬ 
ing a full use of the capabilities of CTA. 

As a summary of our detailed results, charts in Figs. [T71 and 
m compare the performance parameters for the five modes at 
the trigger and analysis levels. 

The parameters are given for low and high energies; for the 
latter £ = 3 TeV is used at both levels. For low energies, we 
choose E around the energy threshold, which is different at both 
levels and then we use E - IQ GeV for the trigger level (Fig. 
[TtI) and E = 150 GeV for the analysis level (Fig.fTSl). 

Dependence on the source position with respect to the axis 
of the array is illustrated by using two values of - 0° and 
4.9°. The former value is less important for the sky scans, be¬ 
cause small 0off contribute negligibly to Acc, however, it al¬ 
lows to illustrate the performance of different modes in deep 
observations on-source. The latter value is crucial for sky scans, 
as the largest contribution to Acc at high energies comes from 


<^off - (4-6)° (except for mode N) and at low energies in modes 
2D and 2C from (^off - (3 - 5)°. For an easier comparison, 
for each value of E and 0oft the best value of a given parame¬ 
ter among the five modes is used to scale the parameters of all 
modes, so that the best value of each parameter corresponds to 
1 . 

As we can see, in both ranges of E the largest collection areas 
at large (^off correspond to mode 2D, which property is reflected 
also in the largest values of Acc for this mode. Mode D has a 
similarly large A(0ofr = 4.9°) at the trigger level, however, it has 
a less efficient gamma/hadron separation at low energies, lead¬ 
ing to a larger loss of gamma-ray events, and then at the analysis 
level its parameters are worse than in 2D. Also at all energies 
2D has a larger collection area at large 0off, ~ 6°, which results 
in larger Acc in this mode. For modes N, C and 2C, collection 
areas at large (pos, as well as acceptance values, are much lower 
which disfavors these modes for sky scans. This conclusion is 
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Figure 17: Summary of performance parameters for the five modes at the trigger 
level for £ = 70 GeV (left column) and 3 TeV (right column). The top (on- 
axis) panels are for (p^fi = 0° and the middle (off-axis) for t^off = 4.9°. For each 
pai'ameter, x, its best value in the five modes, Xb, is used to scale the parameter 
values so the charts show x/jCb for the collection area and acceptance and JCb/x 
for the angular and energy resolution. 


Figure 18: Summary of performance parameters for the five modes at the anal¬ 
ysis level for E = 150 GeV (left column) and 3 TeV (right column). The top 
(on-axis) panels are for t^off = 0° the middle (off-axis) for 0off = 4.9°. 
Values of parameters are scaled by the best value, as in Fig. [m 


further confirmed by the mean sensitivities of scans (Section 
4.2.4), which indicate the best scan performance in modes D 
and 2D and, at intermediate energies, the overall best perfor¬ 
mance of mode 2D. 

In contrast, the parameters for 0off = 0° favor mode C for on- 
source observation. This mode has the best angular and energy 
resolution at high energies, and at low energies its resolutions 
are similar to other modes. In addition, for the energies around 
the analysis energy threshold (a; 150 GeV) mode C offers the 
less biased energy reconstruction (see Fig.fTObl. It also has the 
largest collection area for low E and, although at high E it has 
a small A, the rejection of background events is very efficient in 
this mode for (pos — 0° and then for on-axis observations mode 
C typically gives the best sensitivities. 

The standard mode N appears disfavored for the types of ob¬ 
servation studied in this work. Nevertheless we do not regard 
the N mode as an entirely inefficient. In the wobble mode of 
observation, in which the telescopes are slightly offset from the 
source position, the N mode might appears as an optimal con¬ 
figuration for a particular offsets. Somewhat speculatively we 
note also that for the sky surveys using the full CTA, i.e. includ¬ 


ing also the LST subarray (4 telescopes planned in CTA array 
E) and the small sized telescopes subarray (SST; 32 telescopes 
in array E), a conhguration of each subarray in a different mode 
can be considered. Namely, the 2D mode, offering the most 
efficient detection of sources, has a rather poor reconstruction 
performance, especially at high energies. Then, to allow for 
more accurate spectral and morphological studies of the discov¬ 
ered high-energy sources, the SST subarray could be set in the 
C mode, which at high E has the best reconstruction at any 0off ■ 
An assessment of the feasibility of such a strategy obviously 
requires explicit computations. 
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